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Observational constraints on the energy scale of inflation 
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Determining the energy scale of inflation is crucial to understand the nature of inflation in the 
early Universe. Assuming a power-law power spectrum of primordial curvature perturbations, we 
place observational constraints on the energy scale of the observable part of the inflaton potential by 
combining the 7-year Wilkinson Microwave Anisotropy Probe data with distance measurements from 
the baryon acoustic oscillations in the distribution of galaxies and the Hubble constant measurement. 
Our analysis provides an upper limit on this energy scale, 2.3 x 10 16 GeV at 95% confidence level. 
Moreover, we forecast the sensitivity and constraints achievable by the Planck experiment by per- 
forming Monte Carlo studies on simulated data. Planck could significantly improve the constraints 
on the energy scale of inflation and on the shape of the inflaton potential. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Inflation in the early Universe has become the stan- 
dard model for the generation of cosmological perturba- 
tions in the Universe, the seeds for large-scale structure 
and temperature anisotropies of the Cosmic Microwave 
Background (CMB). The simplest scenario of cosmolog- 
ical inflation is based upon a single, minimally coupled 
scalar field with a fiat potential. Quantum fluctuations 
of this inflaton field give rise to a Gaussian, adiabatic 
and nearly scale-invariant power spectrum of curvature 
perturbations (see Refs. [3, for reviews). This predic- 
tion is strongly supported by CMB observations from the 
Wilkinson Microwave Anisotropy Probe (WMAP) [1,|J]. 

During slow-roll inflation, the potential energy drives 
an exponential expansion of the Universe. Detecting this 
energy scale is crucial to understand how inflation arises 
in a fundamental theory of physics. It is known that the 
amplitude of the power spectrum of gravitational waves 
is directly proportional to the energy scale of inflation Q . 
Hence one could use a determination of the tensor contri- 
bution to the temperature and polarization anisotropics 
of the CMB to determine this energy scale. 

To our knowledge, the first constraint from CMB ob- 
servations on the energy scale of inflation has been ob- 
tained by Liddle @ . More recent constraints from previ- 
ous releases of the WMAP data on the Hubble scale or 
energy scale during inflation have been presented in 0- 
9 1 ]. If one assumes a specific inflationary model, in some 
cases the energy scale can be pinned down, e.g. for a 
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Landau-Ginzburg potential (i"ol |. 

The WMAP collaboration has released the results of 
7- year observations Q. They updated upper limits on 
the tensor-to-scalar ratio, but not on the energy scale of 
inflation. In this work, we place observational constrains 
on the potential energy scale, the first and second deriva- 
tives of the inflaton potential by using the 7- year WMAP 
data with Gaussian priors on the Hubble constant Hq and 
on the distance ratios of the comoving sound horizon to 
the angular diameter distances from the Baryon Acoustic 
Oscillations (BAO). Our analysis assumes that inflation 
is driven by a single, minimally coupled scalar field and 
that fluctuations on the observable scales are generated 
during an epoch of slow roll of this field (the inflaton). 
We obtain upper limits on three potential parameters 
(Taylor coefficients of the potential). It is clear that the 
WMAP data mainly constrain the amplitude of tensor 
modes by the low-/ temperature and polarization data. 
Compared to WMAP, the Planck satellite is designed to 
measure temperature and polarization anisotropies of the 
CMB to higher accuracy We estimate errors of the po- 
tential parameters for the Planck experiment using the 
Monte Carlo simulation approach. As expected, Planck 
could significantly improve the constraints on the energy 
scale of inflation. 

This paper is organized as follows. In Section|Hl adopt- 
ing a Taylor expansion of the inflaton potential, we ex- 
press the power spectra and their spectral indices of 
scalar and tensor modes at the pivot scale fco m terms 
of the values of the potential energy and its first and 
second derivatives at </>o- In Section IlIII we obtain ob- 
servational constraints on the Taylor coefficients of the 
inflaton potential from the 7-year WMAP data. Using 
a Monte Carlo approach, we analyze the sensitivity of 
the Planck experiment w.r.t. these Taylor coefficients in 
Section ITVl Section IV] is devoted to conclusions. 
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II. THE ENERGY SCALE OF INFLATION 

In standard slow-roll inflation, a scalar field <f> slowly 
rolls down its potential V((f>). The condition for infla- 
tion requires that the potential energy of the inflaton 
field dominates over the kinetic energy. A sufficiently 
flat potential for the inflation is required in order to lead 
to a sufficient number of e-folds. On the other hand, 
when scales relevant to current cosmological observations 
cross the Hubble radius during inflation, the change in 
the value of the inflaton field is typically small. Hence 
the inflationary potential can be expanded Taylor 



n a - f = - 
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(1) 



about the point (f> . Here Vo is the energy scale of infla- 
tion, V\ is the potential force term that is balanced by 
the Hubble friction term during slow-roll inflation, and 
V 2 is the effective mass term of the inflaton field at 4>o. 
This approach is different from the so-called reconstruc- 
tion of the inflaton potential [j| , as we do not track the 
evolution of the inflaton; we rather evaluate the poten- 
tial at the pivot point </V Let us stress that we have to 
assume that the Taylor series in Eq. ([1} converges. 

According to the sign of the second derivative of the 
potential, inflationary models can be classified into large- 
field models and small-field models [TTJ] . The former have 
V 2 > while the latter have V 2 < (see Figures [DEI) . 
Linear models with V 2 = live on the boundary between 
large-field and small-field models. However, this classifi- 
cation might be misleading, as the sign of V2 has nothing 
to do with the initial conditions of inflation. The classi- 
fication is inspired by the Mexican hat potential, which 
leads to inflation for two types of initial conditions: false 
vacuum initial conditions lead to V 2 < 0, while the ini- 
tial conditions of chaotic inflation give V2 > 0. But a 
counter example is the recently studied non-minimally 
coupled Higgs field [fJJ] , which in the Einstein frame has 
an effective potential with V2 < 0, despite being a large- 
field model (initial field values beyond the reduced Planck 
mass). Nevertheless, the sign of V2 bears physical con- 
tent. If V 2 > 0, we can continue the potential to the point 
when inflation comes to an end, but if V2 < 0, it is clear 
that some transition has to happen in order to finish the 
quasi-exponetial expansion of the Universe (see Q). 

We assume that the primordial power spectrum is a 
power-law. To leading order in the slow-roll approxima- 
tion, the power spectra of the scalar and tensor perturba- 
tions, and their spectral indices can be written in terms of 
the value of the potential energy and its first and second 
derivatives at 4>q 
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where Mpi is the reduced Planck mass and r is the tensor- 
to-scalar ratio. The amplitude of the power spectra and 
spectral indices are defined at the scale fco, which corre- 
sponds to the value 4>o of the inflaton field when the mode 
exits the horizon during inflation. Note that only the V 2 
term appears in Eqs. ©-©J which means that ±Vi give 
the same inflationary variables (A Sl r, n s , n t ). Hence ob- 
servations such as CMB anisotropics cannot determine 
the sign of the first derivative of the potential. This fact 
can also be understood from the expansion (H]). Actually 
the sign of V\ is unphysical, because changing the sign of 
V\ is equivalent to </> — > — </> and the sign of <f) itself cannot 
be measured. In what follows we will consider the case 
of Vi > 0. 

In terms of the power spectra and the scalar spectral 
index, the coefficients in the Taylor expansion ((TJ) can be 
expressed (to leading order) as (if r > <=f- V\ 7^ 0) 
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which can be used to reconstruct the inflationary poten- 
tial [j| (see Ref. [l| for a review). In this paper, we use it 
to determine the potential parameters when our Universe 
crosses the Hubble radius during inflation. Note that the 
sign of V 2 is determined by the sign of (3r/8 + n s — I), 
which is consistent with the classification of inflationary 
models in the n s -r plane 



III. OBSERVATIONAL CONSTRAINTS FROM 
WMAP 

We assume a spatially flat ACDM model as back- 
ground cosmology and take the power spectrum param- 
eters (A s ,r,n s ,nt) as input parameters. Then the po- 
tential parameters (Vo,Vi,V 2 ) can be derived by the re- 
lation The reason why we choose the power 
spectrum parameters instead of the potential parameters 
as input parameters in Markov chains is that the power 
spectrum parameters are very sensitive to the values of 
the potential parameters especially close to the origin, 
which leads to an extremely slow convergence of Markov 
chains. Analysis is carried out by using the publicly avail- 
able CosmoMC package [TJj], which explores the param- 
eter space by means of Monte Carlo Markov Chains. Be- 
sides the 7-year WMAP data including the low-Z temper- 
ature (2 < I < 32) and polarization (2 < I < 23) data 0], 
we use two main astrophysical priors: the present-day 
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Hubble constant Hq from the magnitude-redshift rela- 
tion of 240 low-z Type la supernovae at z < 0.1 [lij . and 
the angular diameter distances out to z = 0.2 and 0.35, 
measured from the two-degree field galaxy redshift sur- 
vey and the sloan digital sky survey data [15[ . Following 
the WMAP team we choose k = 0.002 Mpc^ 1 as the 
pivot scale of the primordial power spectra. 

Figure [TJ shows the two-dimensional joint marginalized 
constraints on Vq, V\ and V2 at 68% and 95% confidence 
level (C.L.) from the 7-year WMAP data with Gaussian 
priors on Hq and the distance ratios. We find upper limits 
on the potential energy, the first and second derivative of 
the potential: 

Vq 1/4 < 2.3 x 10 16 GcV, Vl lz < 2.7 x 10 15 GcV, 

\V 2 \ 1/2 < 4.5 x 10 13 GeV, (9) 

at 95% C.L.. Thus the upper limit on the energy scale of 
inflation is less than two orders of magnitude lower than 
the reduced Planck mass, which corresponds to the grand 

unified theory scale. We see that V^ 4 > V^ 3 > IV2I 1 / 2 , 
which is consistent with our assumption that the Taylor 
series converges. Moreover, IV2I 1 / 2 is almost five orders 
of magnitudes below the Planck scale, which confirms 
that the potential must be very flat, consistent with the 
assumption of slow-roll inflation. Another observation 
is that there is a strong degeneracy of the upper limits 
of all considered Taylor coefficients on Vq. This stresses 
again how important the identification of the inflationary 
energy scale is for model building. We have also checked 
that Vq, V\ and V2 are not degenerate with any of the 
other cosmological parameters. 



IV. FUTURE CONSTRAINTS FROM PLANCK 

In this section, following the approach described in 
Refs. [Ill [I?} we generate synthetic data for the Planck 
experiment and then perform a systematic analysis on 
the simulated data. First of all, we assume a fiducial 
cosmological model: baryon density ft^h 2 = 0.0227, cold 
dark matter density Vl c h 2 = 0.108, the Hubble parame- 
ter h = 72.4, reionization optical depth r = 0.089, am- 
plitude of scalar perturbation A s = 2.41 x 10~ 9 , scalar 
spectral index n s = 0.961 and tensor-to-scalar ratio 
r = 0.1. Given the fiducial cosmological model, one can 
use a Boltzmann code such as CAMB [l8[ to calculate 
the power spectra for the temperature and polarization 
anisotropics Cf T , C EE , Cf E and C BB . 

We assume that beam uncertainties are small and that 
uncertainties due to foreground removal are smaller than 
statistical errors. For an experiment with multiple chan- 
nels c with different beam width and sensitivity, the noise 
power spectrum N^ x can be approximated as [l9| 




0.2 0.4 



0.8 1 j.; 
V, [10- 3 M p| f 



FIG. 1: Two-dimensional joint marginalized constraints 1 
and 95% C.L.) on the Vb-Vi plane (top panel), the V0-V2 plane 
(middle panel) and the V1-V2 plane (bottom panel) derived 
from the 7-year WMAP data with Gaussian priors on Hq 
and the distance ratios from the BAO. The dashed line corre- 
sponding to linear models marks the border of large-field and 
small-field models. 



E 



exp 



X FWHM 



-1(1 



1)(^fwhm) 



! /(8 In 2) ,(10) 



(c) 

where a x is the root mean square of the instrumental 
noise per pixel for temperature (X = T), E-mode po- 
larization (X = E) and _B-modc polarization (X = B), 
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,(c) 

'FWHM 

sian beam for channel c. Non-diagonal noise terms are 
expected to vanish since the noise contribution from dif- 
ferent maps are uncorrelatcd. For Planck we combine 
only the 100, 143 and 217 GHz HFI channels, with beam 
width #fwhm = (9.6', 7.0', 4.6') in arcminutes, temper- 
ature noise per pixel ctt = (8.2,6.0,13.1) in [iK and 
polarization noise per pixel oe,b = (13.1,11.2,24.5) in 
\iK (see Ref. [2(| for the instrumental specifications of 
Planck) . 

Given the fiducial spectra C* Y and noise spectra 



XY 



SxyN{ , one can generate a random realiza- 



tion of af rn using the following method [T| 
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and are Gaussian- 



distributed random numbers with unit variance. Then 
one can reconstruct the power spectra of the mock data 

by 
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Once simulated data are produced we perform a Monte 
Carlo analysis through the effective % 2 defined as [2l[ 
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(15) 



where / s k y is the sky fraction sampled by the experiment 
after foregrounds removal. For the Planck data we choose 
/sky = 0.65, corresponding to a ±20° Galactic cut. 

We consider the TT, TE and EE power spectra on 
scales with I < 2000. A measurement of the amplitude 
of the primordial gravitational waves would allow a di- 
rect determination of the inflationary energy scale. Both 
the primordial gravitational waves and the weak gravi- 
tational lensing of the -E-mode polarization are sources 
of £>-mode polarization [22|. On large scales, the contri- 
butions to the -B-mode signal mainly come from the pri- 
mordial gravitational waves generated by inflation. The 
effect of the weak gravitational leansing on the B-mode 
power spectrum ultimately dominates on small scales. 
For the fiducial model with r = 0.1 the B-mode signal 
generated by the weak gravitational lensing dominates 
above I ~ 150 and peaks at I ~ 1000. We consider the 



BB power spectrum up to the peak of the spectrum at 
I ~ 1000 since the uncertainty in C BB arising from in- 
strument noise becomes large on smaller scales. 




Planck 
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FIG. 2: Forecast constraints (68% and 95% C.L.) on the Vo-Vi 
plane (top panel), the Vo-Vi plane (middle panel) and the V\- 
V% plane (bottom panel) for the Planck experiment in the case 
of r — 0.1. The dashed line corresponding to linear models 
marks the border of large-field and small-field models. 



Figure [2] shows the forecast constraints on Vo , V\ and 
Vi at 68% and 95% C.L. from the simulated Planck data 
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in the case of 7' = 0.1. As expected, it illustrates dramat- 
ically how Plank can break degeneracies between the po- 
tential parameters if the scalar-to-tensor ratio is tightly 
constrained by Planck. The marginalized la (68%) errors 
on V , Vi and V 2 are 0.099 [10- 2 M P1 ] 4 , 0.18 [10~ 3 M P1 ] 3 
and 0.27 [10~ 5 A/pi] 2 , respectively. As we can see, Planck 
can place strong constraints on the inflationary energy 
scale, which would provide a firm observational link with 
the physics of the early Universe. The primordial gravi- 
tational waves could be detected at > 95% C.L. based on 
our fiducial cosmological model with r = 0.1. In order 
to distinguish the large-field models from the small-field 
models it is important to detect the sign of the second 
derivative of the potential. The analysis of our fiducial 
cosmological model indicates that Planck will not be able 
to distinguish these models. 

What about the forecast constraints from Planck if 
r <^ 0.1? Here we show the case of r = 0.01. The fidu- 
cial values of the other parameters are taken to be the 
WMAP5 maximum likelihood values Figure [3] shows 
the forecast constraints on the potential parameters at 
68% and 95% C.L. from the simulated Planck data in the 
case of r = 0.01. We can see that in this case Planck only 
provides the upper limits on the energy scale of inflation 
as WMAP7 does. It is implied that Planck fails to detect 
the primordial gravitational waves if the tensor-to-scalar 
ratio is less than one percent. As shown in Figure [3J in 
this case "small-field" models would obviously be favored 
by data compared to the case of r = 0.1. It can be un- 
derstood from Eq. @. If r - 0.01 and n s ~ 0.96, the 
second term dominates over the first one in the bracket. 
Therefore, a red tilt of the power spectrum of curvature 
perturbations leads to a negative second derivative of the 
potential. The upper limits achievable by Planck are 

U 1/4 < 1.4 x 10 16 GcV, Vi /3 < 1.1 x 10 15 GcV, 
\V 2 \ 1/2 < 9.8 x 10 12 GeV, (16) 

at 95% C.L.. 



V. CONCLUSIONS 

In this paper we have placed observational constrains 
on the potential energy scale, the first and second deriva- 
tive of the potential by using the 7-year WMAP data, 
combined with the latest distance measurements from 
the baryon acoustic oscillations in the distribution of 
galaxies and measurement of the present-day Hubble con- 
stant from supernova data. A previous upper limit from 
the first WMAP data release, combined with large scale 
structure data from the 2dF galaxy redshift survey found 

U 1/4 < 2.7 x 10 16 GcV at 90% C.L. 0,0. Our new up- 
per limit on the energy scale of inflation is only slightly 
stronger V 1/4 < 2.3 x 10 16 GeV at 95% C.L., and shows 
a degeneracy with the upper limit on the first derivative 
of the inflaton potential, V^ /3 < 2.7 x 10 15 GeV at 95% 
C.L.. Adding information on the power spectrum of large 
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FIG. 3: Forecast constraints (68% and 95% C.L.) on the V -Vi 
plane (top panel), the V0-V2 plane (middle panel) and the Vi- 
V2 plane (bottom panel) for the Planck experiment in the case 
of r = 0.01. The dashed line corresponding to linear models 
marks the border of large-field and small-field models. 



scale structure, the limit would be improved slightly (23} . 
Both upper limits lie at the energy scale of grand unified 
theories and are consistent with scenarios in which infla- 
tion starts very close to the Planck scale. However, a sce- 
nario in which there is just enough inflation [2~H in order 
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to solve the horizon and flatness problems but V <C Mpi 
for all field values cannot be excluded, nor can false vac- 
uum initial conditions be excluded. Latter scenarios can 
be constrained by a lower limit on the reheating temper- 
ature of the Universe. A recent analysis J25| of the 7- year 
WMAP data finds that the so-called reheating parame- 
ter, fully characterizing the inflationary reheating epoch, 
is constrained by data, which yields lower bounds on the 
reheating temperature, typically > 10 2 GeV — 10 2 TcV, 
dependent on inflation scenarios. 

The Planck experiment will soon provide a very ac- 
curate measurements of CMB temperature and polar- 
ization anisotropies. Using the Monte Carlo simulation 
approach, we have presented forecasts for improved con- 
strains from Planck. Our results indicate that the degen- 
eracies between the potential parameters are broken be- 
cause of the improved constraint on the tensor-to-scalar 
ratio from Planck. Besides Planck also EBEX, a balloon- 
borne CMB polarization experiment, has good chances to 
significantly improve our understanding of inflation [26j j . 

In our analysis, we adopt the three-parameter 
parametrization of the potential for the case of a sin- 
gle, minimally coupled inflaton field. To leading or- 
der in the slow-roll approximation, it is consistent to 
expand the potential only to quadratic order, because 
the third derivative corresponds to higher-order slow-roll 
parameters. The advantage of the model- independent 
parametrization is that we could discriminate inflation- 
ary models by detecting the sign of the second derivative 



of the potential. If we drop one of the other assump- 
tions, e.g. by means of a coupling of the inflaton and 
a Gauss-Bonnet term, we change the prediction of the 
tensor-to-scalar ratio [27| and thus the limit on the infla- 
tionary energy scale. 

It is quite surprising that one of the most fundamen- 
tal questions in the context of inflationary cosmology, 
namely what is the energy scale of the observable part 
of the inflationary expansion, has received relatively lit- 
tle attention in the attempt to extract information from 
cosmic data. We hope that we can convince the reader, 
that even without a detection of tensor fluctuations, some 
interesting upper limits on three parameters can be ex- 
tracted and that our analysis turns out to be fully self- 
consistent. 
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